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Abstract: We investigate the short distance behavior of 3 nucleon forces (3NF) de- 
fined through Nambu-Bethe-Salpeter wave functions, using the operator product expan- 
sion(OPE) and calculating anomalous dimensions of 9-quark operators in perturbative 
QCD. As is the case of NN forces previously considered, we show that 3NF have repulsions 
at short distance at 1-loop, which becomes exact in the short distance limit thanks to the 
asymptotic freedom of QCD. Moreover these behaviors are universal in the sense that they 
do not depend on the energy of the NBS wave function for 3 nucleons. 
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1. Introduction 

Realistic nuclear potentials between two nucleons (2N), determined precisely from 2N scat- 
tering data together with the deuteron binding energy, have often been used to study 
nuclear many-body problems. These two-nucleon forces (2NF), however, generally under- 
estimate the experimental binding energies of light nuclei ||l|, |2[ and this fact indicates the 
necessity of taking into account three- nucleon forces (3NF). In addition, a clear indication 
of 3NF is observed in high precision deuteron-proton elastic scattering data at intermediate 
energies [^. 

The 3NF may also play an important role for various phenomena in nuclear physics and 
astrophysics, which include (i) the backward scattering cross sections in nucleus-nucleus 
elastic scattering Q, (ii) the anomaly in the oxygen isotopes near the neutron drip-line Q, 
and (iii) the nuclear equation of state at high density relevant to the physics of neutron 
stars [p. Universal short distance repulsion for three baryons (nucleons and hyperons) is 
also suggested to explain the observed maximum mass of neutron stars Q. 

Despite of its phenomenological importance, a microscopic understanding of 3NF is 
still limited, due to difficulties in studying 3NF experimentally. Pioneered by Fujita and 
Miyazawa the long range part of 3NF has been modeled by two-pion exchange 
which is known to be attractive at long distance. In addition a repulsive component of 
3NF at short distance is introduced in a purely phenomenological way [ p!o|] . 

To go beyond phenomenology, it is most desirable to determine 3NF directly from the 
fundamental degrees of freedom , the quarks and the gluons, on the basis of QCD. Recently 
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the first investigation of this kind has been attempted using lattice QCD simulations, 
where 3NF have been extracted from the Nambu-Bethe-Salpeter (NBS) wave function for 
a specific alignment of 3 nucleons |11, 13 1. The method used there had been previously 



employed to extract nucleon-nucleon potentials (i.e. 2NF)[p!^, p!5| , |T^, 17] as follows. The 
NBS wave function for 2 nucleons is defined by 

ipE{f) = {0\N{x + f,t)N{x,t)\2N,W) , (1.1) 



where |2N, W^) is a QCD eigenstate with energy W = 2ym^ + A;2 with m^v being the 
nucleon mass, E = k'^/m]\i represents the non-relativistic kinetic energy, and is a nucleon 
interpolating operator made of 3 quarks such as N{x) = e"''"^q°-{x)q^{x)q'^{x). The non-local 
but energy independent potential (more precisely the half off-shell T-matrix) is extracted 
from this NBS wave function as 

{E-Ho)^e{^ = J U{f,r')ipE{r')d^r' (1.2) 

where Hq = —V'^/niN- The non-local potential can be expanded in terms of the velocity 
(derivative) with local function as 

U{f,r') = V{f,V)6^{f-r'), (1.3) 

which becomes 

y(f, V) = Voir) + K(r)(Ji • a2 + VT{r)Si2 + VLs{r)L • S + 0{V^) (1.4) 

at the lowest few orders, where r = |r|, o"j represents the Pauli-matrices acting on the spin 
index of the i-th nucleon, S = (ai + 0^2)72 is the total spin, L = f x p is the angular 
momentum, and 

612 = 3 ^ CTl • 0-2 (1.5) 

is the tensor operator. This method has been shown to work well. The central potentials 
at the leading order in the expansion have qualitatively reproduced common features of 
phenomenological 2N potentials: the force is attractive at medium to long distance while 



it has a characteristic repulsive core at short distance. See also refs. [18, 19| for a summary 
of results and recent developments. 

The present authors have investigated short distance behaviors of the 2NF defined 
in the framework mentioned above, using the operator product expansion (OPE) and per- 
turbation theory thanks to asymptotic freedom of QCD pO| , 21, 22] . (See also a similar 
attempt for the solvable models in 2 dimensions ]^3[.) The behavior of the NBS wave 
function ^pE{r) at short distance (r — )• 0) is encoded in the operator product expansion 
(OPE) of the two nucleon operators: 

Nix/2, 0) N{-x/2, 0)^J2 ^^(^) ^'^(O' 0)' (1-6) 
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where {O^} is a set of local color singlet 6-quark operators with two-nucleon quantum 
numbers. Asymptotically the x-dependence and energy dependence of the wave function 
is factorized into 

^e{x) -^Dk{x){0\Ok{0,0)\2N,W) . (1.7) 

k 

Standard renormalization group (RG) analysis gives |21] the leading short distance behavior 
of the OPE coefficient function as 

Dk{x) ^ (in'^Y' dk , (1.8) 

where ly^ is related to the 1-loop coefficient of the anomalous dimension of the operator O^, 
dk is the tree- level contribution of Dj^iJ)) and finally tq is some typical non-perturbative 
QCD scale. Assuming its matrix element does not vanish, the operator with largest RG 



power Vk dominates the wave function (L7) at short distances. We denote the largest 



power by vi and the second largest one by i'2- 

If vi is non-zero, this leads to the leading asymptotics of the s-wave potential of the 
form 

V{r)^- J'. , (1.9) 



r 



(In^) 



which is attractive for u\> Q and repulsive for i^i < 0. 

If vi = 0, the situation is more complicated. The relative sign of the ratio R between 
the leading and the subleading contributions becomes important and we find: 

V{r) w -R ''^ 1 ,^ . (1.10) 

^' r2(lnni)i-^2 ^ ' 

If R is positive, the potential is repulsive, while it is attractive for negative R. A system of 
two nucleons corresponds to this degenerate case. Unfortunately in this case R depends on 



the energy E. In [21] it is argued that in the relevant energy range the relative coefficient 
R is positive, so that the short distance limit of the nucleon potential is repulsive. 

In this paper, we extend the above OPE analysis to the 3NF. The corresponding equal 
time NBS wave function for 3 nucleons is given by 

V'3iv(r,p) = (0|A^(fi,0)iV(f2,0)A^(x3,0)|^3Jv) (1.11) 

where E^^ and \E2,n) denote the energy and the 3N state. We introduce Jacobi coordinates 
R = {xi +X2 + ^3)/3, r = {xi — X2)/2, p = (xs — {xi +X2)/2)/\/3. From this wave function, 
the three nucleon potential is defined by 



-^Vl - + J2 V2N{nj) + VsNFif, P) 

l^r Pp .^j 



V'SAf (^ P) = Emi^3N{r, p) (1.12) 



where V2N{^ij) with rij = Xi — xj denotes 2NF between (i,j)-pair, V3iyF{f,p} the 3NF, 
Pt = Pp = /'^ the reduced masses. 

In sect. |2|, we start with renormalization group considerations and OPE, which are 
relevant for 3NF. The anomalous dimensions of 9-quark operators are computed in sect. |3|. 
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Finally we discuss the short distance behavior of 3NF in sect. For the convenience of 
the reader we give a brief summary of our results here. The OPE analysis shows that the 
3N central potential at short distance behaves as 

I _4/3tree 

V3NF{f,P) ^ 2^ ^ (X (1-13) 

as s = (r)2 + (p)"^ — t- 0, where (3^^'^'^ is given by 

^^r«e = -14/(33 -2iV/), (1.14) 

where Nf is the number of dynamical quarks. Unlike the 2NF where the situation was 
not completely determined by PT alone, it is shown that the 3N potential always has a 
repulsive core. Furthermore it is universal in the sense that it does not depend on details 
of the 3N state used to define the NBS wave function such as the energy of the state. 

2. Renormalization group analysis and operator product expansion for 
3NF 



2.1 Renormalization group equation for composite operators 

In QCD, using dimensional regularization in D = 4 — 2e dimensions, bare local composite 
operators 0^^\x) are renormalized according to 

0^r\^) = ZAB{9,e)0Pix). (2.1) 

Summation of repeated indices is assumed throughout this paper unless indicated other- 
wise. The meaning of the above formula is that we obtain finite results if we insert the 
right hand side into any correlation function of the fundamental gluon and quark fields, 
provided we also renormalize the bare QCD coupling go and the quark and gluon fields. 
For example, in the case of an n-quark correlation function with operator insertion, which 
we denote by Gn\{go,() (suppressing the dependence on the quark momenta and other 
quantum numbers) we have 

g(^^;)(5,^) = Z^B(5,e) V^'(5,6)gl';l,(<70,6). (2.2) 

We recall from renormalization theory that for the analogous n-quark correlation function 
(without any operator insertion) we have 

gi--)(5,M) = Z-"/2(5,6)af (50,6). (2.3) 

The coupling renormalization is given by 

g^o = fi^'Zi{g,e)g'. (2.4) 

The renormalization constant Zi in the minimal subtraction (MS) scheme we are using has 
pure pole terms only: 



Zi{g,e) = l 5— + -V + 0(5), (2.5) 
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where 



1 fll 2_ 1 „ 1 r34, 2 /13 



Similarly for the fermion field renormalization constant, we have 

Zf{9, 6) = 1 - ^ + 0(/) , (2.7) 

XCf - I 

where jpo = „ n with Cp = — and A is the covariant gauge parameter. The gluon 

field renormalization constant is also similar, but we do not need it here. Finally the matrix 
of operator renormalization constants is of the form 

(1) 2 

ZABig, e) = Sab - + 0(/) . (2.8) 

The renormalization group (RG) expresses the simple fact that bare qiiantities are 
independent of the renormalization scale /x. Introducing the RG differential operator 

the RG equation for n-quark correlation functions can be written as 

{p + |7f(<7)} Gt''\9,n) = Q. (2.10) 

Here the RG beta function is 

/5{g) = eg + Poig, e) = eg - 1%^^^^ = -/3og' - /3ig' + 0(/) , (2.11) 

^+2 dg 

where Poig, e) is the beta function in D dimensions and the RG gamma function (for quark 
fields) is 

7Fig) = Mg, e) = 7F0 g' + 0(/) . (2.12) 



It is useful to introduce the RG invariant lambda-parameter A by taking the ansatz 

A = f^ef^3^ (2.13) 

and requiring PA = 0. The solution is the lambda-parameter in the MS scheme (Ams) if 
the arbitrary integration constant is fixed by requiring that for small coupling 

Finally the RG equations for n-quark correlation functions with operator insertion are of 
the form 

{T^ + ^Mg)} <j:^(5,/x)-7AB(5)^i;B^(5,/^)=0, (2.15) 

where 

iAB{g) = -ZAcPoig, 6)^ = Tiiff' + 0(/) . (2.16) 
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2.2 OPE and RG equations 

Let us recall the operator product expansion for p, r — )• 0: 

Oi(f - p/^/3)02(-r - p/V^)0^{2p/^) ~ DB{r,p)OB{^) ■ (2.17) 

We will need it in the special case where the operators Oi,02,03 on the left hand side 
are nucleon operators and the set of operators on the right hand side are local 9- 
quark operators of engineering dimension 27/2 and higher. All operators in ( |2.17| ) are 
renormalized ones, but from now on we suppress the labels As we know, the nucleon 

operators are renormalized diagonally as 

O, = 0(5,e)of , (2.18) 

and we can define the corresponding RG gamma functions by 

7.(5) = Ma.e) ^ = 7f + 0(/) . (2.19) 



For the nucleon operator, 



7!]^ = 24d, d = ^ = 7. (2.20) 

' 32A^^7r2 967r2 ^ ' 



Next we write down the bare version of ( p. 171 ) (in terms of bare operators and bare 
coefficient functions): 

0f\r-p/V2,)0^^\-r-p/V2,)0f{2p/V2,) ^ {r, p) of {^) . (2.21) 

Comparing to ( ^) , we can read off the renormalization of the coefficient functions: 

DB{r,p) = Ci{9,e)Ug,e)Ug,e)Df\r,p)Z^l{g,e) (2.22) 

and using the /x-independence of the bare coefficient functions we can derive the RG equa- 
tions satisfied by the renormalized ones: 

VDB{g,p^.r,p) + DA{g,P,r,p)^AB{g)=^, (2.23) 
where the effective gamma function matrix is defined as 

lAsig) = lAsig) - [71(5) + 72(5) + 73(5)] ^AB ■ (2.24) 

2.3 Perturbative solution of the RG equation and factorization of OPE 

We want to solve the vector partial differential equation ( |2.23| ) and for this purpose it is 
useful to introduce Uab {g) , the solution of the matrix ordinary differential equation 

m ^ UAB{g) = lAcig) UcB{g) (2.25) 
o-g 
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and its matrix inverse UABid)- We will assume that the coefficient has the perturbative 
expansion 

cei+a2=dA 

= E ''"'P'^" [l^^lo"^ +5'^Ar(^"''^'^-^p) + ' (2-26) 

ai+a2=dA 

where = + and tana; = p/r with r = \r\, p = \p\, and Qr,^p are solid angles 
of the vectors r and p, respectively. Here dA = cIa — (^i + d2 + d-^) is the dimension of 
the coefficient function. Note that in the massless theory operators of different dimension 
do not mix. In the full theory quark mass terms are also present, but they correspond to 
higher powers in r and p, and therefore can be neglected. 

We will also assume that the basis of operators has been chosen such that the 1-loop 
mixing matrix is diagonal: 

jABig) = 2/3o (3a 6ab + 0(/) . (2.27) 

In such a basis the solution of (|2.25| ) in perturbation theory takes the form 

UABig) = {Sab + RAB{g)} g'^^^ , (2.28) 

where RAB{g) = 0(51^), with possible multiplicative log (7^ factors, depending on the details 
of the spectrum of 1-loop eigenvalues Pa- 

Having solved ( 2.25| ) we can write down the most general solution of ( |2.23| ): 

D%^''\g- ps, uj, Qr, ^p) = F^^^'iAs, co, Qr, ^p) UAB{g) • (2.29) 

Here the vector F^^"^ is RG-invariant. Introducing the running coupling g as the solution 
of the equation 

fig) = fig) + Hf^s) = In(As) (2.30) 

F^^"^ can be rewritten as 

F^i"^ (As, oj, Qr, n,) = D^i"^ ig; 1, oj, J)p) UABig) ■ (2-31) 
Since, because of asymptotic freedom (AF), for s — )• also ^ — )• as 

F^^"^ can be calculated perturbatively using ( |2.26| ) and ( p. 28 ). 

Putting everything together, we find that the right hand side of the operator product 
expansion ( 2.17] ) can be rewritten: 

Oiir- p/V3)02i-r- p/V3)03i2p/V3) ^ p"^ 

ai+a2=dB 

X F^''^^As,co,Qr,^p)OBiO). (2.33) 
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where 

dB = UBc{9)0c. (2.34) 

There is a factorization of the operator product into pcrturbative and non-perturbative 
quantities: F^^"^(As,a;, O^, Op) is pcrturbative and calculable (for s ^ 0) thanks to AF, 
whereas the matrix elements of Ob are non-perturbative (but s-independent). Note that 
dc = ds- 

2.4 3NF from OPE 

Using results in the previous subsection, the NBS wave function for 3N can be written at 
short distance as 

i^3Nir,p) E r'''p'''D'X''''{9A,oo,nr,np)UAB{gmOBmE3M). (2.35) 

Since a r^ip"^ term produces angular momenta Zi < ai and h < ct2, we can write 

D'X"''{g,l,U;,nr,np) ~ Yl D'X'"''™{9A,^)ya,md^r)Ya,m2i^p), (2.36) 

mim2 

up to less singular terms at short distances. Then the sum of 2N and 3NF potentials 
y2N+3NF is extracted as 

V2N+ZNF{r, p) = y V2N{rij) + V3NF{r, p) = EzN H , ,^ , (2.37) 

tr' friN V3N{r,p) 

where 

- — — r^— -^ = ^2-^ V^-— — 2— -S = cZ2-^ f2 38) 
r"^ dr^ dr ^2 ' p dp^ dp ^ p^' 

Since 

V2r'Yi^(J^,) = V2p^yz^(0p) = 0, (2.39) 

we have 

(V2 + v2)V'3jv(r-',p) ~ ^ ^ r"iy«,^,(O,)p"=^y„2m2(5^p)(0|OB(0)|i?3iv) 

X (d2 + dl) [z^^i-i"2-2(5, 1,u;)C/ab(5)] ■ (2.40) 

Since terms with = dominate in at short distance, contributions from 7^ 
terms to 2N+3NF potentials are suppressed by an r'^^ p^'^ factor, so that they do not 
contribute at short distance. Therefore we consider terms with c?^ = (ai = 02 = 0) 
hereafter and do not write ai dependence in coefficients. We then have 

(V^ + V2)V'3Jv(r,p) ^ 5^(0|Ob(0)|E3jv) x {d^ + 4) [z?a(5, 1,uj)Uab{9)\ ■ (2-41) 

AB 
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In terms of s and u we write 

^^^'+4^ A], (2.42) 



^ ^ ds ds 



duj"^ sin(2a;) doj 



so that the u dependent part gives a contribution unless w = 0, 7r/2, tt, 37r/2, where 
either r = or p = 0. We assume r ^ and p 7^ in our calculation. Since an u 
dependence appears only at 1-loop or higher orders in Da, we can neglect it unless an 
operator Oa which appears at Ia loop has large anomalous dimension such that Pa — 
is larger than other /3b corresponding to operators Ob appearing at tree level. As we will 
see later such operators are absent; it is then enough to consider the tree level contribution 
in Da, so that uj dependent terms in Da can be neglected. The largest eigenvalue among 
operators appearing at tree level is thus denoted by /3a, which corresponds to z^i discussed 
in the introduction for 2N forces. 
We then obtain 

(4 + 4) [DA(9,l,toOf/.4fl(s)] ~ DA:o(4 + 4)g-"^' 

The NBS wave function is dominated by the term with largest Pa- If we assume Pa is 
non-zero, we finally obtain 

IT^N S (— log(As)) 

3. Anomalous dimensions for three nucleons at 1-loop 
3.1 OPE for 3N operators at tree level 

The general form of a gauge invariant local 3-quark operator is given by 

(x) = Bif;ix) = e'^'^q^/{x)q'^'ix)qf{x) , (3.1) 

where a,/?, 7 are spinor, f,g,h are flavor, a,b,c are color indices of the (renormalized) 
quark field q. The color index runs from 1 to A'^c = 3, the spinor index from 1 to 4, and 
the flavor index from 1 to Af. Note that Bj^^^^ is symmetric under any interchange of pairs 
of indices (e.g. -B;^^!^ = -B^q^J) because the quark fields anticommute. For simplicity we 
sometimes use the notation such as F = fgh and F = a/37 as indicated in (3.1). 



The usual nucleon operator which is employed in lattice simulations is constructed 
from the above operators as 

Biix) = (P+4)««, Bif^^{C^,)p,{^r2y' , (3.2) 

where P+4 = (1 + 74)/2 is the projection to the large spinor component, C = 7274 is the 
charge conjugation matrix, and T2 is the Pauli matrix in the fiavor space (for Af = 2) 
given by {iT2y9 = efa. Both C75 and ZT2 are anti-symmetric under the interchange of two 
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indices, so that the nucleon operator has spin 1/2 and isospin 1/2. Although an exphcit 
form of the 7 matrices is unnecessary in principle, we find it convenient to use a (chiral) 
convention given by 



7fc = . n ' = -, n ' = 71727374 = { ^ ^ • (3-3) 






As discussed in the previous section, the OPE at the tree level (generically) dominates 
at short distance. The OPE of 3N operators given above at tree level becomes 

Bl{x + y- z/V3)B^^{x -y- z/V?.)B^{x + 2z/V?.) = B[{x)BI{x)B^{x) + • • • (3.4) 

where + • • • denote higher dimensional operators, which do not contribute at short distance. 
The leading operators couple only to the states with L = (the S-state). 

If we construct local 3N operators at L = from Ba{x) for nucleons (which only involve 
2 different flavors), there is only one such operator, which has 1 = 1/2 and S = 1/2, due 
to the Pauli statistics of nucleons. Explicitly it is given by 

l^trccJa/3a - ^a^p^a^ " ^ a+a,[2,l] + [2,i] ^ ' 

where f ^ g, f,g = u,d and a / /3, a, /3 = 1, 2, 1 = 3, 2 = 4 for the explicit form of the 7 
matrices. Above no summation is taken for / and a. 

3.2 General formula for the divergent part at 1-loop 

As shown in ref. [^], the gauge invariant part of the divergence from diagrams involving 
exchange of a gluon between any pair of quark fields is given by 

, , 1 1— loop,div J p , 1 fg 

q^/{x)q'f{x) =^ To- q''{x)C^q\x) (3.6) 

6 H 



a, 13 



where 

(3.7) 

for either ai,/3i G {1, 2} (right-handed) or ai,/3i G {3, 4} (left-handed), and it vanishes for 
other combinations. 

The operator in eq. ( p.5| ) can be written as a linear combination of simple operators 
[i?i?i?]p^P^P^. According to this 1-loop formula for divergences, such a simple operator 
mixes only with operators [i?i?i?]p^p^p^ which preserve the set of flavor and Dirac indices 
in the chiral basis as 

F1UF2UF3 = FaUFbU Fc, Pi U Pa U = Fa U Tb U Pc- (3.8) 

Note however that such operators are not all linearly independent. In the case of a 2N 
operator, we have the following constraint 

3[55]fifl + 1 = (3-9) 
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which comes from the general identity 

jY^^ai---ajVj,^6i---6jVc 

Here means a simultaneous exchange between the i-th indices of -Fi,ri and the j-th 
indices of i*2 , r2 . This identity can be generalized to 

gOl-'-ajVc £^l---^iVc gCi-'-cjVc _ ^ai---ai_i6jai+i---ajvc 

i,j,k=l 

from which we have 

S[BBB]^lf^f^ - ± [BBB]\^^Pf-^ = (3.12) 

i,j,k=l 

where the i-th index of ABC, the j-th index of DEF and the A;-th index of GHI are 
cyclically interchanged in (ABC, DEF,GHI)[i, j,k]. For example, 

(ri,r2,r3)[i, 1,2] = (/33/3171, 01/3272,030273), 
(ri,r2,r3)[i,2,3] = (73/3171, "2ai72, "3/33/32) • 

Note that the cyclic interchange of indices occurs simultaneously for both Fj and Fi in the 
above formula. Both 2N and 3N identities are incorporated in our calculation. 

3.3 Results of anomalous dimensions for 9 quark operators at 1-loop 

In tables |l|-Q, we give all eigenvalues of the matrix in units of 2d, which were calculated 
and checked independently by Mathematica and Maple programs, for Fi L) F2 L) F3 = 
{f f f f f9999) with f ^ g and all independent combinations of Fi U r2 U Fa. The five digits 
ni, 723, 77,5, ny, ng in the isospin column give the number of representations with isospins 
1/2,3/2,5/2,7/2,9/2, respectively. For example, 32100 means 3 operators with / = 1/2, 
2 with / = 3/2, 1 with / = 5/2 and with / = 7/2, 9/2. 

The results in the tables show some notable patterns. Firstly the eigenvalues 7j/2(i 
are all even integers; this is non-trivial since there appears to be considerable mixing in our 
initial operator bases. Secondly there is a tendency for the operators with larger isospin to 
have smaller (more negative) eigenvalues. Thirdly there are relations between the entries in 
the tables for different Dirac indices; e.g. the isospin degeneracies for the indices 111113344 
and 111122223 in table 1 are identical and furthermore all the corresponding eigenvalues 
related by a common shift of —32. These observations suggest that there is an underlying 
algebraic structure that we have unfortunately not yet been able to reveal. 



j,k=i 



.ai ■■■aj_if)feaj+i---ajV(. iaj6fc+i-- -ftiv, 



(3.10) 
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Note that a combination obtained from an other one by 

i) the interchange (1,2) -H- (3,4) or 

ii) the simultaneous interchange of 1 -f-)- 2 and 3 -f-)- 4 or 

iii) the interchange 1^-7-2 

has obviously the same spectrum of anomalous dimensions and for this reason not listed 
separately. 

The star symbol * next to the eigenvalues means that there is a corresponding operator 
which overlaps with the tree operator in eq. (|3.5D . Since the tree operator has I = 1/2, 
this can only happen if the corresponding ni is different from zero. The tree operator is 
invariant under the symmetry i) above whereas its two spin components are exchanged 
under the symmetry ii). This means that whenever the tree operator overlaps with a 
particular operator, it also overlaps with these transformed ones. The meaning of the 
symbol is an overlap between the tree operator and the transform of the operator under 
symmetry iii). It is intriguing that the tree operator generally overlaps with the largest 
eigenvalue for given Dirac indices; again a fact for which we do not yet have a simple 
explanation. 

As a simple example, let us consider the Dirac index distribution 111112222 (third 
entry in table 1). The part of the tree operator relevant for this case is 

Ti = iB((i - B(if) {Bid - Bil!) {Bfd - Blif ) . (3.13) 

There are altogether 53 local 9-quark operators with this distribution of Dirac indices, but 
the number of independent ones is reduced to 2 after imposing all the gauge identities ( p^ ) 
and ( |3.12| ). A possible choice is 

Oi = b((I bU! Bfll, 02 = Bill Bilf Bill (3.14) 

Using the gauge identities we have 

Ti = - 502, (3.15) 



which is proportional to the operator (of / = 1/2) corresponding to anomalous dimension 
-60 in our units. The other combination 

O1 + 3O2 (3.16) 

has anomalous dimension -84 and 7 = 3/2 and has no overlap with Ti. As explained above, 
among the operators with Dirac index distribution 333334444 there is one with anomalous 
dimension -60 which also overlaps with the spin=l component of the tree operator and 
among the ones corresponding to 222221111 one which overlaps with its spin=2 component. 
The space of operators corresponding to other Dirac index distributions are considerably 
larger than in this example. For example, the case 111223344 (last entry in table 4) has 
1369 operators before, and 117 operators after imposing the gauge identities. 

As can be seen from tables the largest eigenvalue of 7^^^ is 16d (occurring already 
at the tree level). Therefore, the largest eigenvalue of becomes 2(i x (8 — 36), which is 
negative, so that /3^'^*' = —14/(33 — "^Nf). Therefore, in conclusion, the operators at the 
tree level in OPE dominate at short distance in the 3N NBS wave function. 
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4. Short distance repulsion of 3NF 



As discussed before, the 3N potential at short distance is given by 



1 -4/3*?*"^ 
mN s'^[— log(As)j 



where (3^J^^'^ is given in eq. (1-14). Since this result dominates over the one appearing 
in the 2N potential, which is of the form ( |1.1[1| ), the above behavior of V2N+iNF{'f'^ p) 
at short distance must come solely from V^nf{''', p)- Unlike for the 2NF no additional 
nonperturbative considerations are required in this case and therefore we can conclude 
that the 3NF is always repulsive at short distance, which is universal in the sense that it 
does not depend on the details of the 3N state, used to define the NBS wave function, such 
as its energy E. Note however that this conclusion is valid for the 3NF in our definition 
{i.e. defined from the NBS wave function) since potentials are not observable in general 
and therefore scheme-dependent. Unless one fixes the scheme for the definition of the 
potential , it is meaningless to ask whether the 3NF has a repulsive core or not. By using 
the potential scheme considered in this paper, we can show that the 3NF universally have 
repulsive cores. 

Finally we would like to note that as listed in Tables 1-4, the total number of 9-quark 
local operators is several hundred and the spectrum of anomalous dimensions is rather 
dense. We singled out the ones with largest anomalous dimensions which dominate at short 
distance but it is not evident at what distance scale this leading behavior sets in. Even if 
it turns out that these individual operators really dominate at extremely short distances 
only, our main conclusion may remain valid due to the fact that all the eigenvalues of the 
effective gamma matrix are negative (corresponding to short distance repulsion) . We think 
a simple explanation of this fact should exist (maybe related to the Pauli principle). 

An interesting and important extension of the present analysis is to investigate the 
short distance behavior of the three baryon force (3BF) by the same method. Its results 
can tell us whether there is a universal short distance repulsion also in the 3BF, which has 
been suggested to explain the observed maximum mass of neutron stars Q . 
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Table 1: Eigenvalues 7j of the anomalous dimension matrix 7 and isospins of the corresponding eigen- 
vectors for the case 5f4g. The five digits nl, n3, n5, ri7, n9 in the isospins column give the number of 
representations with isospins 1/2,3/2,5/2,7/2,9/2, respectively. 
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